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Abstract. 

& . We investigate the magnetic and glassy transitions of the square-lattice XY 

| model in the presence of random phase shifts. We consider two different random- 

shift distributions: the Gaussian distribution and a slightly different distribution 
I (cosine distribution) which allows the exact determination of the Nishimori line where 

magnetic and overlap correlation functions are equal. We perform Monte Carlo 
simulations for several values of the temperature and of the variance of the disorder 
, distribution, in the paramagnetic phase close to the magnetic and glassy transition 

lines. We find that, along the transition line separating the paramagnetic and the quasi- 
7-H ■ long-range order phases, magnetic correlation functions show a universal Kosterlitz- 

Thouless behavior as in the pure XY model, while overlap correlations show a disorder- 
£^ ■ dependent critical behavior. This behavior is observed up to a multicritical point which, 

, in the cosine model, lies on the Nishimori line. Finally, for large values of the disorder 

C*~) ' variance, we observe a universal zero-temperature glassy critical transition, which is 

— ■ in the same universality class as that occurring in the gauge-glass model. 
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1. Introduction 

The two-dimensional XY model with random phase shifts (RPXY) describes the 
thermodynamic behavior of several disordered systems, such as Josephson junction 
arrays with geometrical disorder [HE], magnetic systems with random Dzyaloshinskii- 
Moriya interactions [3], crystal systems on disordered substrates [I], and vortex glasses 
in high-T c cuprate superconductors [5] . See [HI [7] for recent reviews. The RPXY model 
is defined by the partition function 

Z({A})=eM-n/T), 

U = -yKe^ x U xy il) y = -^Jcos^ - 9 y - A xy ), (1) 

(xy) {xy) 

where ip x = e x , U xy = e tAxy , and the sum runs over the bonds (xy) of a square lattice. 
The phases A xy are uncorrelated quenched random variables with zero average. In most 
studies they are distributed with Gaussian probability 

Pa{A xy ) oc ex P (-^j . (2) 

We denote the RPXY model with distribution ([2]) by GRPXY. We also consider the 
RPXY model with distribution (cosine model) 

P c (A xy ) oc exp (^^j , (3) 

which we denote by CRPXY. Such a model is particularly interesting because the 
distribution fl3]) allows some exact calculations along the so-called Nishimori (N) line 
T = l/P = a In both GRPXY and CRPXY models the pure XY model is 

recovered in the limit a — > 0, while the so-called gauge glass model [TO] with uniformly 
distributed phase shifts is obtained in the limit o — )■ oo. 

The nature of the different phases arising when varying the temperature T and 
the disorder parameter a and the critical behavior at the phase transitions have been 
investigated in many theoretical and experimental works [3j [T], [2], SI El El d El El Ell El 

iniiisiiHEaiiaiiaiiBiiiaE^ 

[351 EHl ET1 [391 EHl HQl SH S21 Sal SH SSI HS1 SZl SU SHI Eol [5U E21 ESI [5H [55]. in spite 

of that, a conclusive picture of the phase diagram and of the critical behaviors has not 
been achieved yet. 

The expected T-a phase diagram for the GRPXY and CRPXY models, which is 
sketched in Fig. [H presents two finite-temperature phases: a paramagnetic phase and a 
low-temperature phase characterized by quasi-long-range order (QLRO) for sufficiently 
small values of a; see, e.g., [55] and references therein. The paramagnetic phase is 
separated from the QLRO phase by a transition line, which starts from the pure XY 
point (denoted by P in Fig. [I]) at (a — 0,T — T X y ~ 0.893) and ends at a zero- 
temperature disorder-induced transition denoted by D at (o~d,T = 0). The QLRO 
phase extends up to a maximum value o~m of the disorder parameter, which is related 
to the point M = (o~m,Tm), where the tangent to the transition line is parallel to the 
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Figure 1. Phase diagram of RPXY models as a function of T and of the disorder- 
distribution variance a. 



T axis. No long-range glassy order can exist at finite temperature for any value of a, 
including the gauge-glass limit a — > oo [2TJ [22] . Several numerical studies of the gauge- 
glass XY model HU E3 E21 SHI SSI H9j [50J [51] support a zero-temperature glassy 
transition. A more complete discussion of the known features of the phase diagram will 
be reported below. 

In this paper we investigate several controversial issues concerning the critical 
behavior at the magnetic and glassy transitions in RPXY models. In particular, we 
will check whether the critical behavior along the paramagnetic-QLRO transition line is 
universal and belongs to the universality class of the pure XY model, whether there is a 
multicritical point along the paramagnetic-QLRO transition line, and, finally, whether 
the T = glassy transition extends from o = oo to a = <td, see Fig. [H and belongs 
to the same universality class as that in the XY gauge-glass model. For this purpose, 
we perform Monte Carlo (MC) simulations of the GRPXY and CRPXY models for 
several values of the temperature and of the variance a, approaching the magnetic and 
glassy transition lines from the paramagnetic phase. As we shall see, our results for 
the CRPXY model provide a robust evidence for a universal Kosterlitz-Thouless (KT) 
behavior of the magnetic correlations along the paramagnetic-QLRO transition line 
from the pure XY point P to the point M where the transition line runs parallel to the 
T axis and magnetic and overlap correlations are equal. Along the line the magnetic 
correlation length £ behaves as ln£ ~ u^ 1 ^ 2 , where u t is the thermal scaling field, and 
the magnetic susceptibility as x ~ £ 7 ^ 4 (corresponding to 77 = 1/4). On the other hand, 
the behavior of the overlap correlations appears to be a dependent along this transition 
line. Moreover, the numerical results for the CRPXY model indicate that the point M 
is multicritical. We conjecture that these conclusions hold for any RPXY model. In all 
cases we expect that the paramagnetic-QLRO transition line is divided into two parts by 
a multicritical point M, where magnetic and overlap correlations have the same critical 
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behavior, though they are not equal. At variance with what happens in the CRPXY 
model, the point M is not expected to coincide with the point in which the tangent to 
the transition line is parallel to the T axis: this coincidence should be a unique feature 
of the CRPXY model. Then, from P to M we expect any RPXY model to behave as the 
CRPXY, that is a KT behavior for magnetic correlations and a a dependent behavior 
for disorder-related quantities. The universality of the behavior has been confirmed by 
our numerical results for the GRPXY model. 

Finally, we have investigated the critical behavior for large values of o. Our 
numerical results provide strong evidence for a universal zero-temperature glassy 
transition for a > cro- For T — > overlap correlation functions are critical, and, in 
particular, the corresponding correlation length £ D diverges as £ G ~ T~ u when T — > 
with v = 2.5(1). 

This paper is organized as follows. In Sec. [2] we review the known results for the 
phase diagram and for the critical behavior of the RPXY models. Sec. [3] provides the 
definitions of the quantities considered in our numerical work. In Sec. H] we study the 
critical behavior along the thermal paramagnetic-QLRO transition line which starts at 
the pure XY point P and ends at multicritical point M. In Sec. [5] we discuss critical 
behavior along the N line of the CRPXY model and show that the point M where the 
N line intersects the critical line is multicritical. In Sec. [6] we investigate the glassy 
critical behavior at T = for o > o~d- Finally, in Sec. [7] we draw our conclusions. There 
are also several appendices. |Appendix A| reports some details of the MC simulations. 



Appendix B is devoted to a careful analysis of the KT renormalization-group (RG) 



equations and of the corresponding RG flow. We derive the most general form of the (3 
function for the sine-Gordon model and discuss the structure of the scaling corrections 
in the XY model. These results are used in the discussion of the behavior at the 



paramagnetic-QLRO transition. In Appendix C we discuss some features of the critical 
behavior at a multicritical point. In Appendix D we briefly discuss the RG equations in 
the presence of randomness. Finally, in Appendix E we report some analytical results 
for the magnetic correlations in the gauge-glass model. 



2. The phase diagram 

In Fig. [H we show the expected T-a phase diagram of the RPXY models. In the absence 
of disorder (a = 0) the model shows a high-T paramagnetic phase and a low-T phase 
characterized by QLRO controlled by a line of Gaussian fixed points, where the spin-spin 
correlation function (^ x ip y ) decays as l/r' 7 ^ for r = \x — y\ — > oo, with rj depending on 
T. The two phases are separated by a Kosterlitz-Thouless (KT) transition [56] at [57] 
Pxy = I/Txy = 1.1199(1). For r = T/T XY - 1 -> 0+ the correlation length and the 
magnetic susceptibility diverge exponentially as ln£ ~ r -1 / 2 and x ~ £ 7 ^ 4 > respectively. 
An interesting question is whether these features change in the presence of random phase 
shifts. 

The low-temperature phase of RPXY models shows QLRO for sufficiently small 
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values of a. The universal features of the long-distance behavior are explained by the 
random spin- wave theory [3], obtained by replacing 

cos(0 x -By- A xy ) ^l--{9 x -9 y + A xy f (4) 

in Hamiltonian (Tj[|). This scenario has been accurately verified by Monte Carlo (MC) 
simulations in both GRPXY and CRPXY models [22]. The QLRO phase disappears 
for large values of a, see, e.g., [6] and references therein; more precisely, as we shall see, 
for cr > 0.31 in the case of the CRPXY model. 

For cr — y oo phases are uniformly distributed and one obtains the gauge-glass 
model. Even if this model has been much investigated [51 [TOj HU [151 EE 113 HI EES 

EH E21 EHl ESI EHl E21 Eal EH ESI EEl STil SH 113 SH SSI SHI IHl SHI SHI EOl EH E21 El], 

its phase diagram and critical behavior are still controversial. No long-range glassy 
order can exist at finite temperature [21} [22] . However, this does not exclude more 
exotic low-temperature glassy phases [ISl HZ], for example a phase characterized by 
glassy QLRO. Many numerical works at finite and zero temperature support a zero- 
temperature transition [51 [332 ESI S31 SSI SSI SSI EDI EI]- According to this scenario, 
the correlation length determined from the overlap correlation function diverges as 
£ D ~ T~ u when approaching the critical point T = 0. The critical exponent v has 
been estimated by finite-temperature Monte Carlo (MC) simulations, obtaining [15] 
\jv — 0.39(3) and [49J 1/V = 0.36(3). The exponent v is related to the T = stiffness 
exponent 9 by 9 = —\jv. The T = numerical calculations of [13] and [51] provided 
the estimates 9 = —0.36(1) and 9 m —0.45 respectively, which are consistent with the 
finite-temperature estimates of v. The T = transition scenario has been questioned 
in [ISl HU HU S3 S3 E21 ESI El], which provide some numerical and experimental 
(for Josephson-j unction arrays with positional disorder [53]) evidence for the existence 
of a finite-temperature transition at T w 0.2, with a low-temperature glassy phase 
characterized by frozen vortices and glassy QLRO. 

Other features of the phase diagram are better discussed within the CRPXY model, 
characterized by the random phase-shift distribution ([3]), because of the existence of 
exact results along the so-called Nishimori (N) line [8j E] 

T = 1/(3 = a. (5) 

Along the N line the energy density E is known exactly: 

where Io{P) and ii(/3) are modified Bessel functions. Moreover, the spin-spin and overlap 
correlation functions are equal: 

[<^y>] = [|<M,)| 2 ]. (7) 

As already noted in [9], the N line should play an important role in the phase diagram, 
because it is expected to mark the crossover between the magnetic-dominated region 
and the disorder-dominated one. 
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Figure 2. Phase diagram of the ± J (Edwards- Anderson) Ising model on the square 
lattice. The phase diagram is symmetric under p — > 1 — p. 



In the GRPXY and CRPXY models, the paramagnetic phase is separated from 
the magnetic QLRO phase by a transition line, which starts from the pure XY point 
(denoted by P in Fig. Q]) at (a — 0, T — Txy ~ 0.893) and ends at a T = transition 
point induced by disorder (denoted by D) at (a D , T = 0), where An important 

result has been proven for the CRPXY model [8] : the critical value ou of a along the N 
line is an upper bound for the values of a where magnetic QLRO can exist. Therefore, 
at the critical point M = (ctm,7m) the tangent to the critical line should be parallel to 
the T axis; moreover, the critical value o\d at T = must satisfy od < om- 

It is worth noting how similar the phase diagrams of the CRPXY model and of the 
square-lattice ±J Ising model in the T-p plane are, see Figs. [I] and [2] respectively. The 
square lattice ± J (Edwards- Anderson) Ising model is defined by the Hamiltonian 

7~t-±J = — JxyO~x°~y, (8) 

(%y) 

where a x = ±1, the sum is over pairs of nearest-neighbor sites of a square lattice, 
and J X y are uncorrelated quenched random variables, taking values ± J with probability 
distribution P(J xy ) = P^(J X y — J) + (1 —p)S(J xy + J)- This model presents an analogous 
N line [58] in the T-p phase diagram, defined by tanh(l/T) —2p+l = 0. The transition 
point along the N line is a multicritical point (MNP) [59J [60] . Moreover, the critical 
behavior for T > Tmnp and T < Tmnp is different. From the pure Ising point at p — 1 

| We mention that the first renormalization-group (RG) analyses based on a Coulomb-gas 
description [3] predicted <jd = 0, but it was later clarified that this was an artefact of the 
approximations. Indeed, experimental [TT] and numerical works [TT1 [50] as well as refinings 
of the RG arguments [H [23j [27j [28l [3_H [35], have shown the absence of a reentrant transition for 
sufficiently small values of a. 
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to the MNP the critical behavior is analogous to that observed in 2D randomly dilute 
Ising (RDI) models [61]. From the MNP to the T = axis the critical behavior belongs 
to a new strong-disorder Ising (SDI) universality class [62]. Finally, the T = end-point 
of the low-temperature paramagnetic-ferromagnetic transition line is the starting point 
of a T = transition line, characterized by a glassy universal critical behavior |63j. 

In [8] it was also argued that, in the RPXY models (in particular, in the CRPXY 
one) the low-temperature paramagnetic-QLRO transition line from the critical point M 
to the point D runs parallel to the T axis, so that o~n — o~m- The same arguments fail in 
the 2D ±J Ising model [621 [59j [60j Ell [65] , although they provide a good approximation. 
Thus, they are likely not exact also in the case of the RPXY models, although they may 
still provide a good approximation, suggesting that < ou — &d °~m- 

In the phase diagram reported in Fig. which refers to the CRPXY, we may 
distinguish two transition lines meeting at point M: the thermal paramagnetic-QLRO 
transition line from P to M, which can be approached by decreasing the temperature 
at fixed a, and the transition line from M to D, which can be instead observed by 
changing disorder at fixed T for sufficiently low temperatures. As we shall see, our 
numerical results for the CRPXY model provide some evidence that the point M is 
multicritical. We conjecture that the same conclusion holds for generic RPXY models, 
though in the generic case we do not expect the multicritical point M to coincide with 
the point where the tangent to the critical line is parallel to the T axis. 

The phase transition from the paramagnetic to the QLRO phase is generally 
expected to be of KT type (ln£ is expected to have a power-law divergence), but 
its specific features, for instance the precise form of the power-law behavior and the 
value of the exponent rj, have not been conclusively determined yet. Some numerical 
results supporting the KT-like behavior were presented in [30J. The disorder-driven 
T = transition at cr D has been argued [231 [2H [301 EHl 02] to show a KT-like behavior 
with ln£ ~ (er — o~d) 1 and x ~ £ 2_r? with rj = 1/16. However, other RG studies 
[31] [28] obtained a different behavior: ln£ ~ (a — cr D ) -1 / 2 ! . The value of r\ associated 
with the magnetic two-point function has been believed to vary along the critical 
line [31 [231 [281 [31], from 77 = 1/4 of the pure XY model at a = to 77 = 1/16 at the T = 
transition. As we shall see, our numerical results along the thermal paramagnetic-QLRO 
transition line, from P to and including M, strongly support 77 = 1/4, independently of 
a. 

In the following sections we investigate some of the open issues of the RPXY 
models, by performing MC simulations of the GRPXY and CRPXY models close to their 
magnetic and glassy transition lines. In particular, we investigate the critical behavior 
at the thermal paramagnetic-QLRO transition line (from point P to the multicritical 
point), along the N line in the CRPXY model, and at the T = glassy transition line 
for large disorder. 
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3. Notations 

We consider RPXY models defined on square lattices of size L 2 with periodic boundary 
conditions. We define the magnetic spin-spin correlation function 

G(x -y) = [(V»x W)] (9) 
and the overlap correlation function 

Go (x-y) = [\(^ y )\ 2 }. (10) 

The angular and square brackets indicate the thermal average and the quenched average 
over disorder, respectively. The latter can also be written in terms of the overlap 
variables. Given two copies of the system with spins ipx and ipi 2 \ we define 

g x = ^Vi 2) , G (x-y) = [(q x q y )], (11) 

where the thermal average is performed over the two systems with the same disorder 
configuration. We define the magnetic susceptibility x = J2 X G(%), the overlap 
susceptibility Xo = J2 x Go(x), and the second-moment correlation lengths 

^ 2 ^ G(0) - G(g min ) ^ 2 _ G o (0) -G (g min ) 

^min^^rnin) (?min^° (Omin) 

where g min = (2ir/L,0), q = 2 sing/2. 
We also define the quartic couplings 

^~2W' ^^^[(H 4 )-2(|//| 2 ) 2 ], (13) 

922 = X22 = p: (K|^| 2 ) 2 ] - [(H 2 )] 2 ) , (14) 

9c = 94 + 3fi-22, (15) 

where fi = J^x^x and V = L 2 . Note that for the pure XY model g 2 2 = and g c = g±. 
Finally, we define an overlap quartic coupling g as 

9o = X4o = ^[(l/iol 4 )] - 2[(|/. | 2 )] 2 , (16) 

where \i a = Y, x Qx- 

4. Critical behavior along the thermal para-QLRO transition line 

In this section we study the critical behavior of the RPXY models along the thermal 
paramagnetic-QLRO transition line, see Fig. (TJ which starts at the point P on the 
a = axis and ends at the multicritical point, which belongs to the N line in the 
CRPXY model. For this purpose, we perform MC simulations of the GRPXY and of 
the CRPXY model for several values of T and a in the paramagnetic phase, where the 
magnetic correlation length £ is large but finite. Fig. [3] shows the points where the 



simulations are performed. The MC algorithm is described in Appendix A We average 
over a large number of samples, A^ s 10 4 in most cases. We consider large lattice sizes, 
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Figure 3. Values of T = 1//3 and a where MC data were collected. The circles 
and crosses refer to the GRPXY and CRPXY models, respectively. The dotted line 
T = a is the N line for the CRPXY model. We also show some estimates of T c for 
the GRPXY and CRPXY models, and the critical point (MNP) of the CRPXY model 
along the N line. The dashed line is the prediction (|?7|) for the behavior of T c at small 
values of a. 



satisfying L/C, ?Z 10, in order to make finite-size effects negligible and obtain infinite- 
volume results. The residual finite-size effects are in all cases smaller than, or at most 
comparable with, the statistical errors. 

In the following we first discuss the critical behavior of the magnetic spin-spin 
correlation function fl9]). We show that disorder is apparently irrelevant: for any o 
the correlation length diverges following the KT law valid for o = and the magnetic 
susceptibility diverges with critical exponent r\ equal to 1/4. Then, we discuss the 
behavior of observables related to the overlap correlation function ([10]) . finding that the 
critical behavior of these quantities is apparently a dependent. 

4-1- Critical behavior approaching the pure XY transition point 

We wish now to understand the critical behavior along any line that lies in the 
paramagnetic phase and ends at the pure XY critical point at o = and T = Txy- For 
a = 0, as T approaches the critical temperature T X y from above (paramagnetic phase), 
the magnetic correlation length £ diverges as 

ln(£/X) = Cr- 1 ' 2 + 0{r 1 ' 2 ), t = (T - T XY )/T XY , (17) 

where X and C are nonuniversal constants. In the case of the square-lattice XY model 
with nearest-neighbor interactions [57] (3xy = l/^xr = 1.1199(1), X = 0.233(3) and 
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Figure 4. MC estimates of £ for /3 — Pxy = 1-1199 and several values of a versus 
CT -i/2 rpj^g hashed ij ne corresponds to a linear fit to ln£ = C CT cr -1 / 2 + b. 
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1.776(4)3 The magnetic susceptibility \ diverges as, see |Appendix B 

h 



X 



Htm 



+ O(l/ln 2 



Note that while A x is a nonuniversal amplitude, the coefficient b x of the leading 



logarithmic corrections is universal. As shown in Appendix B, it can be computed 
from the perturbative expansion of the RG dimension of the spin variable, obtaining 

b x = 7T 2 /16. 

We now consider the GRPXY model and study the critical behavior of x an d £ as 
one approaches the pure XY critical point along the line (3 = (3xy = 1.1199 by decreasing 
a. We collected data for 0.46 ^ a <; 0.14 in the infinite- volume limit, corresponding to 
the quite large range of correlation lengths 4 < £ < 50. Fig. 0] shows a plot of In £ versus 
(j-i/2 The data fall on a straight line, showing that for a — > 

\n£~o-- 1/2 . (19) 

This behavior can be understood within the RG framework. The general discussion 
presented in Appendix C shows that, as long as disorder is less relevant than the thermal 
perturbation, the critical behavior can be simply obtained by replacing r with the 
nonlinear thermal scaling field. Note that it is not necessary that disorder is irrelevant 
to obtain the result ( fl9l) . In general, the thermal nonlinear scaling field Ut is an analytic 
function of the system parameters. Thus, in the presence of disorder it is a function of 

§ Equation (TIT)) holds whatever the definition of the correlation length is, but of course X depends on 
the specific choice for £. Reference [57] studied the exponential correlation length £ gap , which is defined 
as the inverse of the mass gap, and determined the corresponding constant X gap = 0.233(3). Since 
in the critical limit [55] £ 2 /£g ap = r = 0.9985(5), the constant X for the second-moment correlation 
length we use is given by X = X gapV /r = 0.233(3). 
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both r = (T — Txy) /Txy and a such that, close to the XY transition point, it behaves 
as 

u t (T,a) = t + c a a + . . . (20) 

where the dots stand for higher-order terms. If disorder is less relevant than the thermal 
perturbation, then 

Ht/x) = Cu; 1/2 + o(ul /2 ), (21) 

along any straight line in the T, a plane which ends at the XY pure transition point. 
Since this relation also holds for a = and Ut(r, 0) = r, C and X are the same constants 
reported below ( TiTl) . Along the line T = Txy Equation ( 12TI) implies 

Hi/X) = -^- TJ - 2 +0{a 1 l% (22) 
\c a o~ j 

in agreement with the observed behavior. In order to determine c a we have performed 
fits to 

ln(£/X) = C>- 1/2 (l + H, (23) 

using X = 0.233(3). We obtain the estimates C a = 2.010(2) and b w -0.11. In 
particular, a fit of the data satisfying £ > 7 gives C a = 2.0102(8) and b = —0.108(2), 
with x 2 /DOF w 1.1 (DOF is the number of degrees of freedom of the fit). Using 
C = 1.776(4) and C a = C/^/c^, we obtain 

= 0.781(4). (24) 

The constant c a is nonuniversal and as such is model dependent. However, for o — > 
the fields A xy are typically very small and the distribution functions for the GRPXY 
and CRPXY models are identical to leading order in A xy . We thus expect that the first 
correction to the thermal scaling field due to disorder is identical in the two models, i.e. 

u t ,GKPxy(r, a) = u t ,cKPXY(r, o) + 0(cr 2 ), (25) 
which implies that c a is the same in the GRPXY and CRPXY models. 



4-2. Critical behavior of the magnetic correlations at fixed a 

Standard arguments that apply to critical lines and multicritical points imply that the 
critical temperature at fixed a must be the solution of the equation 

u t [T c (a),a} = 0. (26) 

Therefore, Equation ( 1201) also implies that for small values of a the critical temperature 
for the GRPXY model (and also for the CRPXY model if (125]) holds) is given by 

T c (a)=Tx Y {l-c a a + 0(a 2 )]. (27) 

Equation (1271) can be checked by analyzing data at fixed small values of a. We have 
performed MC simulations of the GRPXY model at a — 0.0576 for several values of {3, 
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Figure 5. Plots of ln£ vs t" 1 / 2 , where t = (T — T c )/T c , for the GRPXY and CRPXY 
models at a — 0.0576. For both models we use T c — 0.8528, as obtained by using (|2"7)l. 
The dashed line corresponds to a fit of the GRPXY data to ln£ = dr 1 !" 1 + a. The 
dotted line that connects the MC data is drawn to guide the eye. 



from f3 = 0.95 to ft = 1.02, corresponding to 10 < £ < 26, and of the CRPXY model 
at the same value of a for f3 = 0.92, 0.95, 0.99 corresponding to 7 < £ < 16. In Fig. [5] 
we plot £ versus t~ 1 ^ 2 with t = T/T c — 1 and T c = 0.8528 given by (j27p [if we take the 
errors on T X y and c CT into account, we have T c = 0.8528(3)]. Clearly, £ — > oo as t — > 0, 
confirming ( 1271) . Moreover, they are clearly consistent with the KT behavior 

\n£ = ar 1/2 + b. (28) 

A fit of all available data for the GRPXY model to (EHJ) gives a = 1.841(2) and 
b = -1.511(5) (with x 2 /DOF w 1.3) keeping T c = 0.8528 fixed. A nonlinear fit, taking 
T c as a free parameter, gives T c = 0.852(2), in good agreement with ( 127|) . Note that the 
estimate of the constant b is close to the corresponding XY- model value lnX = —1.46(1). 
This is no unexpected since X(a) = X + 0(a). 

We also collected data at a = 0.1521 for both the GRPXY and CRPXY models, for 
0.8 < f3 < 1.1199 (corresponding to 2 < £ < 37) and 0.96 < (3 < 1.145 (corresponding 
to 5 ^ ( ^ 46), respectively. Again, the data fit well the KT behavior ( ]2~8j) . see 
Fig. E Fits of the MC data for £ > 10 to (EED (for which x7 D F < 1) give the 
estimates T c = 0.772(2) for the GRPXY model, and T c = 0.762(1) for the CRPXY 
model. Note that ( 1271) would give T c = 0.7872 for o = 0.1521, which is slightly larger 
than the above estimates. This is not unexpected since, when increasing a, higher-order 
terms (which are different for the two models) may become important in ( 1201 . We also 
mention the estimates b = -1.82(7) and b = -1.78(3) for the GRPXY and CRPXY 
model, respectively, from which one obtains estimates of the corresponding length scale 
X(a) = e b , X = 0.162(11) and X = 0.169(5). We also determined £ for other values 
of a, but in a smaller range. The results are compatible with a KT behavior, but they 
do not allow us to get robust estimates of T c . We only mention that in the case of the 
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Figure 7. Plot of ln(x/£ 7/4 ) versus 1/ln^/X. We fix X = 0.233, which is the length- 
scale value valid for the pure XY model. We show data for the pure XY model (taken 
from [67]), and for the GRPXY and CRPXY models at various values of cr, at T = Txy 
and along the N line. The dashed line corresponds to a fit to a + 7r 2 /(16 \n£/X) of the 
pure-XY data satisfying £ > 10 (we obtain a = 0.8058(1) with x 2 /DOF « 0.7). 



GRPXY at a = 0.1936, for which we have only data for £ < 20, we find T c w 0.74. 

At a KT transition the magnetic susceptibility behaves as in f fl8|) . where b x = 7r 2 /16 
is universal. In Fig. [7] we show x/£ 7 ^ 4 f° r the GRPXY and CRPXY and several values 
of o together with those of the pure XY model taken from [ST]. We report the data 
versus ln£/X(cr = 0). We could have also used ln£/X(cr), where X[o) is determined 
from the fit of £. This choice gives a plot essentially identical to the one reported, which 
is not unexpected since, by using ln£/X(cr = 0) or ln£/Y(cr) one simply changes the 
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estimate g* = g\ = 13.65(6) obtained by form-factor calculations [57] . 



corrections of order a/ In 2 which are present anyway. The results appear to follow 
the same curve within the errors (except those obtained along the N line, which we shall 
discuss in Sec. [5]). They provide strong evidence that the value r\ = 1/4 is universal 
along the thermal paramagnetic-QLRO transition line. Also the slope appears universal 
(the coefficient b x does not depend on a), as expected on the basis of the discussion of 
Appendix B. The constant A x corresponds to the intercept of x/£ 7 ^ 4 a ^ m £/A(cr) = 0. 



As it can be seen from the figure, this constant, which is not universal, varies very little 
with a: differences are not visible within our errors, except for the CRPXY data at 
a = 0.307. However, note that for this value of a the critical behavior is controlled by 
the multicritical Nishimori point, i.e. by the special point M which appears in Fig. [TJ 
we will return to it in Sec. [5j 

In conclusion, the above numerical results provide a strong evidence that the 
magnetic two-point correlations show a KT behavior along the thermal paramagnetic- 
QLRO transition line in GRPXY and CRPXY models. 

4-3. Quartic couplings 

We now discuss the behavior of the quartic couplings defined in (fl3j) - (fT5l) . We recall 
that in the pure XY model #22 = while g^ = g c behaves as 

^4 = ^+ (ln ^ x)2 +Q(l/ln 4 ; (29) 



where gl and b g are universal; see Appendix B We mention the estimates g\ = 13.65(6) 



obtained by form-factor computations in [57], and g\ = 13.7(2) by field-theoretical 
methods [68]; other results for g\ can be found in [69J and references therein. 
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Figure 9. Estimates of <?22 versus l/£ for the GRPXY model at fixed j3 = /3xy = 
1.1199. The line is a fit of g 22 to c£ _1 . 

In Fig.Elwe show some MC results of g c for the CRPXY model at a = 0.1521, 0.0576 
and the GRPXY model at (3 = (3xy — 1-1199 (within our errors of a few per mille the 
infinite-volume limit is reached for L/£ ?t 10, as in the pure XY model [67]), and 
compare them with MC results for the pure XY model taken from [67]. The results are 
identical within errors. For example, if we consider the CRPXY model for a = 0.1521, 
a fit to g* + b g /(\n£/X) 2 gives g* = 13.57(10) and b g = -3.1(1.4), with x 2 /DOF « 0.4, 
to be compared with the value [67] g\ = 13.65(6) of the pure XY model. Both g* and 
b g , which are universal in the pure-XY universality class, do not depend on a. 

The quartic coupling g 2 2 defined in f ll4|) is interesting because it is particularly 
sensitive to randomness effects, since in the pure XY model it vanishes trivially. The 
estimates of (722 in the GRPXY model for T = Txy and several values of a are shown 
in Fig. [9j They decrease with decreasing a, and appear to vanish when a — > as 

222 ~ cr e , (30) 

with e w 1.0. A fit to fl30]) gives e = 0.97(4), c = 3.1(3) with x 2 /DOF w 1.1, where 
DOF is the number of degrees of freedom of the fit. 

The fast decrease of #22 along the line T = Txy [note that g 2 2 ~ l/£ implies 
#22 ~ exp(— ccr -1 / 2 )] might suggest the irrelevance of disorder, and therefore that the 
critical value g% 2 vanishes along the thermal paramagnetic-QLRO transition line. This 
conclusion is apparently contradicted by the results at fixed a > 0. The results for the 
CRPXY model at various values of a, a = 0.0576, 0.1521, 0.2992, 0.307, are shown in 
Fig. [TOl where they are plotted versus (ln£/X)~ 2 , which is the correction expected in 
the pure XY model for RG invariant quantities. The coupling g 2 2 is quite small, but 
definitely different from zero on the transition line. For a = 0.1521 an extrapolation 
using + 6/(ln£/X) 2 suggests a nonzero critical limit. Using only data satisfying 
f > 10, this fit gives g* 22 = -0.068(8) and b = -0.080(15), with x 2 /DOF w 0.4. We 
should also mention that the data for the largest values of £, those satisfying £ > 10 say, 
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show linear extrapolations to the critical point. The data denoted by a plus along the 
line related to a — 0.0576 is obtained by using ([50")) and ([25]) with the results of the 
fits along the T = Txy line. 



may be consistent with a vanishing critical limit, but only assuming a slower logarithmic 
approach, i.e., #22 ~ b/(ha.£/X). For instance, the data with £ > 10 are consistent with 
this behavior (the fit gives b = -0.482(4) with xV D0F « 1-1). At a = 0.0576 the 
l/(ln£) 2 extrapolation of the data satisfying 7 ^ £ ^ 16 gives g% 2 = —0.008(6) with 
X 2 /DOF w 1.3. The data of g 2 2 at a ~ 0.30 are larger, but this can be explained by 
crossover effects, since this value of a is quite close to the critical point along the N line, 
where the critical behavior may change, see Sec. [5j 

Overall the results for g 2 2 suggest a nonuniversal critical value. 



4-4- Critical behavior of the overlap correlations 

We now discuss the critical behavior of overlap correlations, cf. (II II) . which are the 
appropriate quantities to understand the role of disorder. We consider the critical 
behavior of the overlap susceptibility which is expected to behave as Xo ~ £,o~ Vo - m 
the case of the pure XY model we have t] = 2t] = 1/2. In [55] it was noted that the 
following relations 

2rj - r] w - for GRPXY, (31) 

7T 

1 2 

2rj - r] ss - — ^— for CRPXY (32) 

7T 

approximately hold in the whole QLRO phase (within the small statistical errors), even 
very close to the KT transition, as long as a is not to large (in practice a should not 
be close to o~m, where M is the Nishimori point defined in Fig. CD). This would suggest 
that they may remain valid up to the transition. Given the strong numerical evidence 
that the exponent 77 associated with the magnetic correlation is rj — 1/4, see Sec. 14.21 
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the above relations imply that rj varies along the paramagnetic-QLRO transition line 
approximately as 

r] « - - - for GRPXY, (33) 

2 7T 

_ 1 a + a 2 /2 for CRpxY; ( 34 ) 

2 7T 

at least for sufficiently small values of a . We wish now to verify if the high-temperature 
data are consistent with these predictions. In Fig. [11] we plot Xo/£, 2 ~ Vo versus 1/ ln(£/X). 
The scaling is reasonable. We also report Xo/^ 2 ' 71 ", fixing rj to the pure XY value 
r] = 1/2. Again the ratio is consistent with a limiting finite value. However, if 
Xo behaves as in the pure XY model, we would expect a a-independent slope, see 



Appendix B which is not supported by the data. 



We now consider the ratio £ /£ between the second-moment correlation lengths 
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obtained from the overlap and spin correlation functions, cf. (I12p ]||l In order to estimate 
this ratio in the case of the pure XY model, we performed MC simulations (using the 
cluster algorithm) in the range 0.93 < (5 < 1.033 corresponding to 12 < £ < 110. Taking 
into account the logarithmic scaling corrections, i.e. fitting the XY-model data satisfying 
£ > 32 to a + b/(\n£/X) 2 with X = 0.233, we obtain the estimate £ /£ = 0.417(4). In 
Fig. [12] we show the results for several values of a. They are all consistent with a finite 
critical value, confirming that the paramagnetic-QLRO transitions are characterized 
by a single diverging length. The results can be extrapolated by assuming = 
a + b/(\n£/X) 2 for £ ->■ oo. We obtain £ /£ = 0.417(5), 0.428(5), 0.425(7), 0.425(3) for 
the GRPXY model at a = 0.0576, 0.1521, 0.1936 and the CRPXY model at a = 0.1521, 
respectively. A larger result is found for the CRPXY model at a w 0.299, 0.307: 
£o/£ « 0.49. 

These results indicate that the ratio varies along the transition line, although it 
changes very weakly for small values of a. Again, this is consistent with the observation 
that disorder- related quantities, like i] and #22, depend on a. 



5. Critical behavior along the N line in the CRPXY model 

We now consider the critical behavior along the N line T = o in the CRPXY model, 
approaching the transition point from the paramagnetic phase. We recall that along 
the N line the magnetic and overlap correlation functions are equal, so that t] = rj and 
£ = £ exactly. We performed several MC simulations along the N line, in the range 
1.5 < /3 < 2.4, corresponding to 2 < £ < 28, and considered large lattice sizes, in order 
to obtain infinite- volume results. 

Our MC estimates of the magnetic correlation length £ are consistent with an 

|| In a Gaussian theory without disorder, in which the magnetic correlation function is given by 
G(p) = (p 2 + m 2 )- 1 , one can easily find that = = 0.408248... 
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exponential increase, i.e. with a behavior of the form ln£ ~ t 1 I 2 with t = T/Tm — 1, 
see Fig. [6j A linear fit to 



with x 2 /DOF < 1.0. We also mention that alternative fits to £ = at~ b and to 
ln£ = at -1 + b give rise to significantly larger x 2 ■ 

In order to estimate the exponent 77, we fit x an d £ to x = c £ 2_r? - Considering 
the MC results satisfying £ > £ min = 5, we find rj = 0.246(4) with xVDOF ~ 1-0. If 
we increase £ min , r\ slightly decreases, but it is always compatible with r] — 1/4. These 
results suggest that r\ = 1/4 also along the N line. 

Fig. [13] shows the estimates of g c . The critical limit of g c is consistent with the 
results for the pure XY model and those obtained along the thermal paramagnetic- 
QLRO line at smaller values of a, see Fig. El Indeed, a fit of all data of g c to ( 1291) gives 
g* = 13.49(13) with x 2 /DOF w 0.6. If we consider only the data satisfying £ > 4, we 
obtain g* = 13.6(3). 

The above-reported results (KT behavior of £, rj = 1/4, and g* ~ glxy) suggest 
that the magnetic correlations behave as in the pure XY model. There is, however, a 
result which contradicts this hypothesis. As we discussed in Sec. 14.11 the rate of approach 
of x£~ 7//4 to its limiting value, should be universal. As can be seen from Fig. this is not 
the case: the slope of the data along the N line is clearly different from that predicted 
for the pure XY model. Thus, even though at the Nishimori point the magnetic critical 
behavior is the same as that observed along the thermal paramagnetic-QLRO transition 
line, corrections are different, implying the presence of a new (probably marginal) RG 
operator, which only contributes to scaling corrections in magnetic quantities. 

A better evidence for the presence of a new, disorder-related operator is obtained 
by considering g 22 and £ /£. In Fig. [13] we also report estimates of g 22 along the N line 
and along the line a = 0.307. If the estimate ([36]) is correct, the two lines intersect the 
critical line at the same point, the Nishimori point. It is quite clear from the data that 
the limiting value of #22 along the two lines is quite different. A fit of all available data 
on the Nishimori line to g* 22 + &/(ln£/X) 2 gives g* 22 = -7.00(5) with xV D0F ~ °- 9 - 
On the other hand, a fit of the data along the line at fixed a = 0.307 gives g 22 ~ —0.8. 
The same phenomenon is observed for the ratio £ D /£- As can be seen in Fig. U21 for 
a = 0.307 this ratio is approximately equal to 0.49, which is clearly different from the 
result that holds exactly along the Nishimori line, £ G /£ = 1. The large differences of 
the values of these two RG invariant quantities along the two lines provide compelling 
evidence that the Nishimori point is a multicritical point as in the 2D ±J Ising model 

To understand this conclusion, let us review the basic results that apply to 
multicritical points. The singular part of the free energy should obey a scaling law 



ln£ = at~ 1/2 + b 
of the data satisfying £ > 5 gives the estimate 
T M = o-m = 0.307(2), 



(35) 



(36) 



•FsingCui,^) =b d T sing (b yi Ui,b y2 u 2 ) 



(37) 
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Figure 13. Estimates of g c and <?22 along the N line and at a = 0.307. The dotted 
lines indicate the estimate [67] g* = g\ XY = 13.65(6) for the pure XY model. 



where u\ and u 2 are two relevant scaling fields. They can be inferred by using the 
following facts: (i) the transition line at M must be parallel to the T axis, since it has 
been proved [8] that % is an upper bound for the values of a where QLRO can exist; 
(ii) the condition T = a at the N line is RG invariant. We therefore have 

U\ — a — um + ■■■ (38) 

where the dots indicate nonlinear corrections, which are quadratic in Act = a — o~m and 
AT = T — Tm-i so that the line u% — runs parallel to the T axis at M. Moreover, we 
choose 

u 2 = T - a, (39) 

so that the N line corresponds to u 2 = 0. 

Close to the multicritical point, any RG invariant quantity, such as #22, is expected 
to behave as 

R = f R {u lU - y ^). (40) 

Now, the N-line corresponds to u 2 = 0, so that a RG invariant quantity converges to 
/r(oo). On the other hand, the line o = oyi corresponds to u\ = 0, so that a RG 
invariant quantity converges to /^(0) which is generically expected to be different from 
(00). Thus, if the Nishimori point is multicritical, we expect RG invariant quantities 
to have a different critical value along the two lines. This is exactly what we observe 
for g 22 and £ /£. Thus, in view of the numerical results we conclude that the Nishimori 
point is a multicritical point. 

It is interesting to note that the multicritical behavior is not observed in the 
magnetic sector. For instance, g* along the paramagnetic-QLRO line is equal to its 
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XY value g\ X v The same result holds along the Nishimori line. In terms of the scaling 
function f 9c defined in f l4"0]) these results imply 

/*(0) = /*(oo) = gtxr • (41) 

It is then natural to conjecture that g* = g\ XY along an y hne that intersects the 
Nishimori point, i.e. that f 9c {x) = g\ XY for any x. The absence of multicritical behavior 
in the magnetic sector is also supported by the fact that £ always shows a KT behavior 
and that the magnetic exponent rj at the Nishimori point is equal to the pure-XY value 
1/4. 

The results we have presented should apply to generic RPXY model. In all 
cases we expect a multicritical point along the paramagnetic-QLRO transition line. It 
follows from universality that, at the multicritical point, the magnetic and the overlap 
correlation functions have the same critical behavior — hence, we have rj = rj — though 
they may not be necessarily equal as is the case for the CRPXY model. Note that 
this point is not expected in general to coincide with that in which the tangent to the 
transition line is parallel to the T axis. 

6. Glassy critical behavior at T = 

In the limit o — > oo the RPXY model corresponds to the gauge-glass model in which the 
phase shifts are uniformly distributed. This model has been extensively studied both 
at zero and at finite temperature [SI [TTil [HI [IHl [TBI [13 [IHl [ISl Ell 

[MlEgEEllISiDllalS^ [2U [22] showed that no 

long-range glassy order can exist at finite temperature. Although this result does not 
exclude the possibility of a finite-temperature transition with an exotic low-temperature 
glassy phase, for example a phase characterized by glassy QLRO, most numerical works 
PH Uni EEJ H31 05J H6j HHJ EOJ [5Tj support a zero-temperature glassy critical behavior. The 
overlap correlation length £ diverges as T~ v for T — > 0. We mention the estimates [45] 
1/u = 0.39(3) and [49J \/u = 0.36(3) from finite-temperature MC simulations, and [43] 
\jv = 0.36(1) and [51J 1/z/ ~ 0.45 from T = numerical calculations. Moreover, if one 
assumes that the ground state is nondegenerate in the overlap variables, one obtains 
that at T = the finite-size overlap susceptibility satisfies the relation Xo = L 2 , so that 
r] = 0. We mention that this scenario was questioned in [4Til4^l4lll47ll4^[5^1,[53l [54]. 
which claimed the existence of a finite-temperature transition at T ^ 0.2. 

A natural scenario for the phase digram of the GRPXY and CRPXY models is that 
the glassy transition, which occurs for a = oo, is not isolated but that it is the endpoint 
of a phase transition line that starts at the paramagnetic-QLRO transition line. In 
particular, if the zero-temperature glassy transition scenario applies to the gauge-glass 
model, we expect a line of T = glassy transitions for any a > ctd, see Fig. [TJ A natural 
conjecture would be that all these transitions belong to the same universality class. 

To check this scenario we performed MC simulations of the CRPXY model at 
o~ = 2/3, 5/9, 1/2, oo, which are larger than op < a M w 0.31. As we shall see, the 
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Figure 14. MC estimates of the correlation lengths £ and £ for the CRPXY model 
at a = 2/3. 



results clearly support a glassy T 
of the gauge-glass model. 



transition in the same universality glass as that 



6.1. MC simulations 

We performed MC simulations of the CRPXY model on square L x L lattices with 
periodic boundary conditions. Most of the results we shall present refer to runs with 
o~ = 2/3. In this case we considered L = 20, 30, 40, 60, 80 and temperatures between 
T = 2/3 (at the Nishimori line) and T = 0.1 (for L = 80 we considered 0.22 < T < 2/3). 
We averaged over a relatively large number N s of samples: N s = 6000, 9000, 7000, 
3000, and 2000 samples for L = 20, 30, 40, 60 and 80, respectively. We used the 



MC algorithm discussed in Appendix A combined with the parallel-tempering method 
[70| ITT] . Moreover, to check the universality of the transitions, we also performed 
parallel-tempering MC simulations for a = 5/9 and lattice sizes L = 60, 70 (5000 and 
1000 disorder samples, respectively), a = 1/2 and L = 70 (1000 samples), and o = oo 
and L = 20, 30, 40, 60 (5000, 5000, 2000, 2000 samples, respectively). The points in 
the T-a plane where we collected MC data are shown in Fig. [3j 

At the glassy transition the critical modes are those related to the overlap variables, 
while the magnetic ones are noncritical. This is clearly shown in Fig. [T3J which shows £ 
and £ D for a = 2/3. The overlap correlation length £ D increases steadily with decreasing 
the temperature, while the magnetic correlation length £ freezes at sufficiently low 
temperatures at a value £ ~ 3.3. Therefore, the critical temperature and exponents 
must be determined from quantities related to the overlap correlation functions. 
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Figure 15. MC estimates of the ratio R$ a = £ Q /L for the CRPXY model at a = 2/3. 
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Table 1. Estimates of T c obtained by fitting to g2]) with n = 6. DOF is the 
number of degrees of freedom of the fit. 



6.2. Evidence for aT = glassy transition at a = 2/3 

In order to determine the critical temperature, we analyze R^ o = £ /L. The results, 
shown in Fig. [151 show no evidence of a crossing point in the range of values of T of the 
data, T > 0.1, and thus provide the bound T c < 0.1 for the critical temperature T c . A 
more precise determination of T c can be obtained by a finite-size scaling (FSS) analysis. 
We fit the data to 

%=P„[(T-T c )L 1 /-] > (42) 

keeping T c and v as free parameters. Here P n (x) is a polynomial in x of order n. The 
order n is fixed by looking at the x 2 °f the fit. For each n we determine the goodness 
X 2 (n) of the fit. Then, we fix n such that x 2 { n ) is n °t significantly different from 
X 2 {n + 1)- The results we report correspond to n = 6. To identify the role of the 
corrections to scaling we repeat the fit several times. Each time we fix two parameters 
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Table 2. Estimates of v obtained by fitting R^ o to (gSJ) with T c = and n = 6. DOF 
is the number of degrees of freedom of the fit. 

T max and L min and we only include the data which correspond to lattices satisfying the 
conditions T < T max and L > L m[n . 

In Table [H we report the estimates of T c for several values of T max and L min . We 
obtain estimates of T c which are quite small and satisfy the upper bound 

T c < 0.01 . (43) 

Since our data satisfy T > 0.1, this estimate allows us to conclude that our results are 
fully consistent with a zero-temperature transition. From now on, we always assume 
T c = 0. 

6. 3. The critical exponent v 

In order to determine the critical exponent v related to the divergence of the correlation 
length £ , we repeat the fit (jl2il at a = 2/3 setting T c = 0. The results are reported in 
Table [2j They slightly increase as T max or L min is lowered, but these changes are small 
compared to the statistical errors. 

In fit (142p we made two approximations. First, we neglected the nonanalytic scaling 
corrections, which decrease as . The results indicate that these corrections are small: 
at fixed T max < 0.25 the estimates of v obtained setting L min = 30 and L min = 40 
differ by much less than the statistical errors. Second, we approximated the thermal 
nonlinear scaling field ut by ut ~ T, neglecting the analytic corrections (see [72] for an 
extensive discussion of this type of corrections). To understand their quantitative role, 
we performed fits to 

R 6 = PniurL 1 ^), u T = T + P T\ (44) 
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Table 3. Estimates of v and p obtained by fitting R^ o to (jUJ) with n — 6. DOF is the 
number of degrees of freedom of the fit. 

where p is a new free parameter. The results are reported in Table [3j Corrections are 
tiny and we estimate \p\ ^ 0.2, so that \ut — T\/T is at most 0.10, 0.02 for T = 0.5, 
0.1, respectively. The estimates of v do not vary significantly and, for L > 30 and 
^max < 0.2, are fully consistent with those obtained before. We quote 

v = 2.5(1) , l/v = 0.40(2) (45) 

as our final estimate. 

To show the quality of our FSS results in Fig. [16] we plot versus TL 1 ^, using 
the estimate (j4"5"l) . All data fall on top of each other with remarkable precision. 
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Table 4. Estimates of r\. On the left we report the results of the fits to (|4"8")l with n = 6, 
on the right those to (|49p with n — 6 and m = 2. In both cases we fix v = 2.5(1). The 
reported errors are the sum of the statistical error and of the variation of the estimate 
of r\ as v changes by one error bar. DOF is the number of degrees of freedom of the 
fit. 

6.4- The critical exponent T) 

As discussed at length in [72], the overlap susceptibility behaves in the critical limit as 

Xo = ulL 2 -v°f( UT L 1 /») . (46) 

Here ut is the temperature nonlinear scaling field, while u h is related to the external 
overlap-magnetic scaling field Uh associated with the overlap variables by Uh = huh(T) + 
0(h 2 ). We have already checked that the thermal scaling field ut can be effectively 
approximated by ut = T. Thus, neglecting nonanalytic scaling corrections, the data 
should behave as 

\n X o = (2- Vo )\nL + \nu h (T) 2 + \nf(TL 1 / u ) . (47) 

We now estimate rj from the analysis of the data at a = 2/3. In a first set of fits we 
set Uh = 1 and approximate In /(re) with a polynomial in x of order n, i.e., we perform 
fits to 

ln Xo = (2 -v) InL + P n {TL l '»). (48) 

The analysis of the x 2 °f the fits indicate that n = 6 allows us to describe accurately 
the data. We fix v to the estimate ( 1451) to avoid an additional nonlinear parameter in 
the fit. The results are reported in Table HI We observe a significant change of the 
estimates as T max decreases; moreover, the quality of the fit is quite poor. This can be 
explained by the presence of sizeable analytic corrections, which means that Uh is poorly 
approximated by a Uh = 1 in our range of temperatures. The same phenomenon occurs 




in the three-dimensional Ising spin glass [72], where the analytic corrections cannot be 
neglected in the analysis of the overlap susceptibility We thus perform a second set of 
fits in which we take into account the magnetic nonlinear scaling field. If we approximate 
lnu\ with a polynomial of order m, we end up with the fitting form 

ln Xo = (2 -v) InL + P n (TL^) + Q rn {T) , (49) 

where we assume Q m (0) = 0. In the following we take m = 2 and again fix v to 
the estimate (j4"5l) . The results are reported in Table HI The quality of the fit is now 
significantly better, indicating that the analytic corrections are important. The scaling 
function is reported in Fig. [TTJand indeed it varies significantly in the range of values 
of T we are considering. The estimates of r\ do not show any systematic variation with 
T max and are always consistent, within errors, with r] = 0. Quantitatively, our data 
allow us to set the upper bound 

\rj \ < 0.05. (50) 
6. 5. Results for the gauge-glass model 

In order to check universality we also performed runs at a = oo, although in this case 
we considered smaller lattices and the errors are significantly larger (partly because of 
the smaller number of samples, partly because of larger sample-to-sample fluctuations). 
The data were analyzed as we did in the a — 2/3 case. First, we determined the critical 
temperature T c . A fit of £ /L to ( H2|) gives rather small estimates of T c . For L min = 20 
we obtain T c = 0.030(2) [0.020(3)] for T max = 0.4 (resp. 0.3). Thus, we can conclude 
that T c < 0.02, which is clearly consistent with T c = 0, given that our data belong to 
the range T > 0.1. The claim that T c w 0.2 is not consistent with our MC data. 
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Figure 18. MC estimates of g a vs T at a = 2/3 for £ = 40, 60, 80. The dotted lines 
correspond to the infinite- volume critical (T = 0) estimate g* = 13.0(5). 



Then, we determined v by assuming T c = 0. The results of the fits to (|42p show 
a significant dependence on T max . For L min = 20, v varies between 2.50(1) and 2.80(4) 
as T max varies between 0.4 and 0.16. If analytic scaling corrections are included, i.e. 
we fit the data to (J4"4"|) . we observe a significantly smaller dependence on T max , but, on 
the other hand, a rather large dependence on L mm , with rapidly increasing error bars 
as L min increases. This is probably due to the fact that we have a somewhat large 
statistical error on the results with the largest value of L, L = 60. The estimates of v 
vary between 2.8 and 3.7 if we take L min = 20, 30 and 0.2 < T max < 0.5 and thus give 
the final result v = 3.3(5). This result is somewhat larger than the estimate (|45p . but 
certainly not inconsistent. It supports — very weakly, though — universality. A better 
check is presented below. 

6. 6. The quartic coupling g Q and universality 

We computed the overlap quartic coupling g Q defined in (fTB"]) . MC results at o = 2/3 are 
shown in Fig. [181 The infinite- volume limit, within our statistical accuracy, is apparently 
reached when L/£ ^ 7, corresponding to T > 0.3 for our largest lattices L = 60, 80. 
The infinite-volume results are quite stable with respect to T, so that we can reliably 
estimate the critical (T = 0) value g*. We obtain 

9*o = 13.0(5). (51) 

According to standard RG arguments, g Q has a universal FSS limit as a function of 
R^ o = £ /L, that is 

g (T,L)=f(R €o ), (52) 

where the function f(x) is universal and satisfies /(0) = g*. This scaling behavior is 
nicely supported by the data at a = 2/3 for various lattice sizes, see Fig. [191 Universality 
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Figure 19. g vs R^ o = £„/L: data at cr = 2/3 for various lattice sizes 
L = 20,40,60,80 (above), and including (below) also data for other values of a: 
a = 5/9, 1/2, oo. 



can be checked by also considering the results for a = 5/9, a = 1/2 and a = oo. Clearly, 
all points fall on top of each other. Note that here there are no free parameters to fiddle 
with and thus this comparison provides strong support to the hypothesis that all these 
models belong to the same universality class. Given the very good evidence we have 
that the model with cr = 2/3 undergoes a T = glassy transitions, this result further 
confirms (and provides stronger evidence than that given in the previous paragraph) 
that the gauge-glass model does not have a finite-temperature exotic glassy transition. 



6. 7. Behavior of the magnetic correlation functions 

Let us now consider the magnetic quantities. The magnetic correlation length £ is zero 
in the gauge-glass model, see Appendix E and increases as one approaches the QLRO 
region. In particular, at T = 0.159, which is below the critical temperature Tm ~ 0.31 
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along the Nishimori line, we obtain f = 3.3(1), 6.7(4), 9.8(3) at a = 2/3, 5/9, 1/2, 
respectively. They are roughly consistent with a behavior like ln£ ~ (a — cr c )~ K assuming 
cr c pa ou ~ 0.30, i.e., with a KT-like behavior along the transition line that connects the 
Nishimori critical point M, see Fig. [TJ and the T = transition point at o = o~d, which 
is expected to run almost parallel to the T axis. Note, however, that while our data 
suggest a power-law divergence of ln£ (therefore, £ has an exponential divergence), they 
are not sufficiently precise to allow us to estimate the power k. The KT value k — 1/2 
is consistent with the data, but k = 1 would be equally reasonable. 

It is also interesting to discuss the behavior of the quartic couplings g c , g 4 , and g 22 



defined from the magnetic correlation functions in ( uBl - ljToj) . In | Appendix E[ assuming 
universality, we predict that, in the critical limit, g± and g 22 should diverge as while 
gdo 2 should go to zero. 

Numerical estimates of g c are shown in Fig. |20j The results are clearly consistent 
with a finite T = limit. Note that the estimates obtained for o = 2/3, 5/9, and 1/2 
are close to the XY value g\ X Y = 13.65(6); actually, they are consistent within errors, 
even at small T, below Tm ~ 0.31. These results are suggestive of a KT behavior of the 
magnetic correlation functions also along the disorder paramagnetic-QLRO transition 
line from M to D, see Fig. [TJ Indeed, for a = 2/3, 5/9, 1/2 we have £ « 3, 7, 10, 
so that along these lines one should be able to observe the critical behavior that arises 
when one approaches the paramagnetic-QLRO transition line at a point with T < Tm- 
In other words, these results imply that the critical limit of g c (o~,T) at fixed T < T M 
along the paramagnetic-QLRO transition line is consistent with the KT value. This 
fact provides some evidence that also along the disorder-driven transition line magnetic 
correlation functions behave as in the pure XY model. Of course, as a increases (thus, 
the magnetic correlation length £ decreases), g c changes significantly and, for a = oo, 
g c is infinite for any T and L. 

The couplings g 22 and g^ are instead expected to diverge as In Fig. [21] we 
report g 22 for the different models. The data are clearly diverging as £ — > oo, but the 
asymptotic behavior g 22 ~ ^ is not clearly observed, likely because the values of £ Q are 
not sufficiently large. Indeed, we only observe that g 22 behaves as Q with k rapidly 
increasing with £ G . More precisely, if we only include data satisfying £ Q < 10 we obtain 
k pa 1. If instead we fit the data with 10 ^ £ G ^ 20 (we have infinite- volume data only 
up to £ D fa 20) we obtain k ~ 1.5. 



7. Conclusions 



We have studied the magnetic and glassy transitions of the square-lattice XY model in 
the presence of random phase shifts and, in particular, the GRPXY and CRPXY model 
defined by the distributions ([2} and d2J). The latter is very useful because it allows 
some exact calculations along the Nishimori line T = a [5J [9], where, in particular, the 
magnetic and overlap two-point functions are equal. We present MC for the GRPXY and 
CRPXY models for several values of the temperature and of the parameter a controlling 
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Figure 20. MC estimates of the quartic couplings g c . The dotted line corresponds 
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the disorder, approaching the magnetic and glassy transition lines from the paramagnetic 
phase. We substantially confirm the phase diagram shown in Fig. [TJ 
Our main results are the following. 

(i) We have carefully investigated the critical behavior along the transition line 
separating the paramagnetic and QLRO phases, from the pure XY point P to the 
multicritical point, which, in the CRPXY model, lies on the N line and is such that 
the transition line runs parallel to the T axis. The magnetic observables show a o~- 

— 1/2 

independent KT behavior: the magnetic correlation length behaves as ln£ ~ u t , 
where u t is the thermal scaling field, u t ~ T — T c (a), and the magnetic susceptibility 
as x ~ £ 7//4 (corresponding to r\ — 1/4). Moreover, the quartic coupling g c defined 
in ( I15p appears to be universal. We obtain g* w 13.6, which is nicely consistent 
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with the corresponding value g\ XY = 13.65(6) of the pure XY model [SZl ESI- 
We have also verified the universality of the leading logarithmic correction to the 
critical behavior of x- On the other hand, the critical behavior of disorder-related 
quantities, such as those related to the overlap correlation function, depends on a. 

(ii) In the CRPXY model, the Nishimori point M, see Fig. [H is a multicritical point 
which divides the paramagnetic-QLRO line into two parts: a thermally-driven 
transition line (from P to M) and a disorder-driven transition line (from M to 
D). This result should be general: a multicritical point should also exist in generic 
RPXY models, although in this case it is not expected to coincide with that 
where the transition line runs parallel to the T axis. Such a multicritical point 
is characterized by the fact that, at criticality, magnetic and overlap functions have 
the same critical behavior, that is rj = r] : in the CRPXY model the two correlation 
functions are exactly equal (more generally, they are equal on the whole N line), but 
we do not expect this property to be generic. It is interesting to observe that the 
multicritical behavior is only observed in the disorder-related quantities. Magnetic 
observables behave, as far as the leading behavior is concerned, as in the pure XY 
model: the correlation length shows a KT behavior, t] = 1/4, and g* = g\ X y m the 
whole neighborhood of the multicritical point. However, corrections are different 
from those appearing in the pure XY model, providing additional evidence for the 
presence of an additional (probably marginal) RG operator, which is responsible 
for the multicritical behavior. 

(iii) Little is known about the behavior along the transition line from the multicritical 
point to D. However, the fact that purely magnetic observables behave as in 
the pure XY model both along the thermally-driven transition line and at the 
multicritical point make us conjecture that the magnetic behavior is also unchanged. 
We have presented some very weak evidence in Sec. 16.71 

(iv) We have investigated the critical behavior for large values of a. We find no evidence 
of a finite-temperature transition for all values of a we have investigated: the system 
is paramagnetic up to T = 0, where a glassy transition occurs. Morever, in all cases 
we verify universality. We can thus conjecture that the critical behavior along the 
whole line that starts in D, see Fig. (TJ is universal: for any a > o~d, one has the 
same critical behavior characterized by the exponents: 

u = 2.5(1), l/u = 0.40(2), \rj \ < 0.05. (53) 

Our estimate of u is consistent with earlier estimates obtained by MC simulations of 
the gauge-glass XY model, for examples l/u = 0.39(3) and l/u = 0.36(3) obtained 
in [45] and [49] respectively, and by numerical calculations of the stiffness exponent 
at T = 0, for example l/u = 0.36(1) and l/u ~ 0.45 obtained in [45] and [51] . Our 
result for i] is consistent with a general argument which predicts rj = 0. 
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Appendix A. Details on the Monte Carlo simulation 

In the simulation we use both Metropolis and microcanonical local updates. The latter 
do not change the energy of the configuration and are defined as follows. Consider a 
site i; the corresponding field is ipi- The terms of the Hamiltonian that depend on ipi 
can be written as 

U t = Re $ iZ ) , Z = Y1 ' ( A - X ) 

3 

where the sum is over all nearest neighbors j of site i. Then, define 

^ = 2^Re(^z)-^ (A.2) 

One can verify that |^| = 1 and that 

Re(fe)=Re(fe). (A.3) 

Thus, the update ipi — > ipl does not change the energy and can therefore be always 
accepted. This update does not suffer the limitations of the Metropolis update: ipi and 
tp'i are not close to each other. 

In our simulation a MC step consists of 5 microcanonical sweeps over all the lattice 
followed by one Metropolis sweep. For each disorder sample we typically perform O(10 5 ) 
MC steps. In some simulations of the CRPXY model we also use the parallel tempering 
method [70l [TTJ . It allows us to obtain results for small values of T, in particular below 
the Nishimori line T = a. In the parallel-tempering simulations we consider Nt systems 
at the same value of a and at N T different inverse temperatures (3 min = {3\, . . . , /3 max , 
where /3 max corresponds to the minimum value of the temperature we are interested in. 
The value /3 min is chosen so that thermalization at (3 = /3 min is sufficiently fast, while the 
intermediate values are chosen so that the acceptance probability of the temperature 
exchange is at least 5%. Moreover, we require that, for some i, Pi = a. This allows us 
to collect data on the Nishimori line. The exact results valid on it allow us to check the 
correctness of the MC code and perform a (weak) test of thermalization. Thermalization 
is checked by verifying that the averages of the observables are independent of the 
number of MC steps for each disorder realization. 

The overlap correlations and the corresponding \o and £ D are measured by 
performing two independent runs for each disorder sample. Finally, note that the 
determination of #22 defined in (I14j) requires the computation of the disorder average 
of products of thermal expectations. This should be done with care in order to avoid 
any bias due to the finite length of the run for each disorder realization. We use the 
essentially unbiased estimators discussed in [731 H2] • 

Appendix B. The KT RG equations 

In this Appendix we consider the RG flow for the sine Gordon (SG) model, with 
the purpose of understanding its universal features. As a results we shall obtain the 



Magnetic and glassy transitions in the square-lattice XY model with random phase shifts3A 



critical behavior of the correlation length and of the magnetic susceptibility at the KT 
transition. This appendix generalizes the results presented in [TUITJECE]- The SG model 
is parametrized by two couplings, a and 5 — we use the notations of [TjJ EHj — whose 
functions are 

a = 2a5 + -^-a 3 + (B.l) 
64 

R s = —a 2 - — a 2 5 + ..., (B.2) 
^ 32 16 ' y J 

where the dots indicate higher-order terms. To all orders the j3 functions have the 

generic form 

a = 2a8 + Kn m a n 5 m , (B.3) 

n+m>2 

Ps = ~k a2+ S hnma n 5 m . (B.4) 

n+m>2 

In the SG model the sign of a is irrelevant, which implies the symmetry relations 

a (a, 5) = -/3 a (-a, 8), Ps{a,5) = s (-a,6). (B.5) 

As a consequence, b a ^ nm = if n is even and bs } nm = if n is odd. Moreover, for a = 
the theory is free and 5 does not flow. Hence 

0s{a = 0, 5) = 0, (B.6) 

which implies bs, nm — if n — 0. 

Let us now consider a general nonlinear analytic redefinition of the couplings 

a = a a , 10 u+ 2j aa,nmM n ^ m , (B.7) 

n+m>2 

5 = asM v + as,nmU n v m . (B.8) 

n+m>2 

We have verified up to the 7 th order that with a proper choice of the coefficients a a ^ nm 
and as^nm one can rewrite the (3 functions in the form 

n (u,v) = - uv, (B.9) 
P v (u, v) = - u 2 (l + b x v + b 3 v 3 + b 5 v 5 + ...). (B.10) 

The couplings u and v are not uniquely defined and indeed there is a family of 
transformations that do not change the functions flB.9j) and fIB.lOj) . Extending the 
previous results to all orders, in the following we assume that we can choose u and v in 
such a way that u (u, v) is given by ( IB. 91) and v (u, v) has the form 

P v (u,v) = -u 2 [l + vf(v 2 )], (B.ll) 

where f(v 2 ) is an analytic function in the region v < vq, where vq is the starting point 
of the RG flow, and satisfies 1 + vf(v 2 ) > in this domain (if this were not true, we 
would have another nontrivial fixed point). This parametrization is unique (universal) 
in the sense that there is no analytic redefinition of the couplings which allows one to 
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write the /3 functions in the form (1B.9|) . (IB. lip with a different function f(v 2 ), i.e. with 
different coefficients b 2n +i- The perturbative calculations of [7j] allow us to determine 

bi = -\. (B.12) 

The analysis of the flow in the general case is analogous to that presented in [7H [76] . 
First, we define the RG invariant function 

Q(u,v) = u 2 -F(v), (B.13) 

F ^ = 2 / n T< 2^ =v 2 + v 3 + -v 4 + 0(v 5 ), 
Jo l + wf{w 2 ) 8 

which satisfies 

where / is the flow parameter. The RG flow follows the lines Q = constant. It is thus 
natural to parametrize the RG flow in terms of Q and v(l). Since 
du 

- = P v (u, v) = -IQ + F(v)][l + vf(v% (B.15) 

we obtain 

l = ~L [Q + F(w)}[i + wf(w 2 )y (R16) 

where v(l = 0) = v . 

Let us now apply these results to the XY model. Repeating the discussion of [771178] 
the XY model can be mapped onto a line in the (u, v) plane with v > 0. The KT 
transition is the intersection of this line with the line Q = and the high-temperature 
phase corresponds to Q > 0. Thus, Q plays the role of thermal nonlinear scaling field, 
i.e. 

Q = qiT + q 2 r 2 + ... (B.17) 

where r = (T — T X y) /Txy- 

To derive the expected critical behavior we consider the singular part of the free 
energy in a box of size L. It satisfies the scaling equation [79] 

J- Sing (r, L) = e~ 2l f(Q, v{l), e~ l L) , (B.18) 

where we have parametrized the flow in terms of Q and v(l) and we have neglected 
all irrelevant operators. If Q > 0, as discussed in [71], v(l) decreases continuously and 
v(l) — > — oo as I — > oo. Since v , the starting point of the flow, is positive, we can fix I 
be requiring 

v(l) = -1, (B.19) 

so that 

r° dw 

1 = / \n^v( Mho. T( 2M = 1 (Q> ( B - 2 °) 
J -i [Q + F(w)][l + wf(w 2 )] 
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It follows 

J- Sing (r, L) = e- 2I ^f(Q, -1, e'^^L) , (B.21) 

which gives the scaling behavior of the free energy (using Q ~ r). In the scaling limit the 
finite-size dependence can be parametrized in terms of where £ is the correlation 
length. This allows us to identify 

£( T ) = &e /(Q ' Vo) , (B.22) 

where £o is a constant. The behavior of £(r) for r — >■ is obtained by expanding I(Q, vo) 
for Q — > 0. The generic behavior is 

/(Q, v ) = -j= W + Yl W(^o)Q n - (B.23) 

The nonanalytic terms in the expansion depend only of the coefficients b 2n +i which 
appear in ( IB.lOp . The first two coefficients are 

h = 7T, 

/! = — = —. (B.24) 
4 16 V 7 

Correspondingly, we obtain 

£(r) = Xexp(7r/^/g)[l + J 1V /Q + O(Q)]. (B.25) 

Expanding Q in powers of r we obtain the celebrated KT expression for the correlation 
length. 

Let us now consider the behavior of the susceptibility. Perturbation theory gives 
for the scaling dimension of the spin correlation function [71] 

1 =-~ + ~5-~5 2 + h 1 a 2 + ..., (B.26) 

where hi is an unknown coefficient. If we perform the redefinitions (a, 5) — > (u, v) 
considered before, we can rewrite 7 a^jjl 

j=----v-— v 2 + ... (B.27) 
' 4 8 16 V ; 

without the a 2 term. In the infinite-volume limit the susceptibility satisfies the scaling 

law 

" rv ® 7(10) + 1/4 



X r 7/4 = ^exp 



dw 

v 



G x [Q,v(l)}, (B.28) 

Pv 

the integral is computed at fixed Q with /3 V given by (IB.15j) . and G x is an analytical 
function. Setting v(l) = —1 and expanding the integral in powers of Q, we obtain an 
expansion of the form 

x t-V4 = A(l + c iy /Q + C2Q + ...). (B.29) 

% The possibility of cancelling the term of order a 2 is related to the existence of a family of 
transformations transformations, given at second order by v! = u + Auv, v' = v + Au 2 with arbitrary 
A, which leave invariant the j3- functions (|B.9j) and (jB.lOj) . By properly choosing A one can eliminate 
the a 2 term in j(u',v'). 
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The coefficient c\ can be computed exactly using the perturbative results (IB. lOj) . (IB. 12j) . 
and (1B.27|) . obtaining 

71 

Cl= 16- 
Using ( IB. 251) we can write 

71 



(B.30) 



ln£/X 



+ 0(ln-^) 



and obtain 



-7/4 



1 + 



7T 



+ 0(1/ In 2 



(B.31) 



(B.32) 



161n(£pO 

Note that the leading logarithmic scaling correction has a universal coefficient. We 
should note that in [76J it was incorrectly claimed that C\ = and, as a consequence, that 
the leading scaling corrections in ( IB. 321) are proportional to l/(ln£) 2 . We numerically 
checked ( 1B.32|) by fitting the infinite-volume numerical data of [67] (more precisely their 



data for (3 > 0.92, corresponding to 10 < £ < 420) to 



ln(xr 7/4 ) = a + 



(B.33) 



obtaining a = 0.804(2) and b = 0.627(9) (with xV D0F ~ °- 7 )> wnich is perfectly 
consistent with the value of b obtained in perturbation theory, i.e. b = rr 2 /16 ~ 0.617 
(fixing b = 7T 2 /16, we obtain a = 0.8058(1) with x 2 /DOF w 0.7, while a fit to 
a + b/ln{£/X) + c/ln 2 (^/X) gives a = 0.8046(9), c = 0.029(22) with x7 D F ~ 0.6, 
which confirms that the next-to- leading correction is very small in (IB. 331) ). 

The result (IB~32l is general. If O is a generic long-distance quantity which behaves 
as in the critical limit, we expect 0/C to behave as x/£ 7//4 > i- e - to satisfy a relation 
analogous to (1B.28j) . It is only needed to replace j(u, v) + 1/4 with the appropriate 
subtracted scaling dimension. Thus, 0/££ also has an expansion of the form ( IB. 321) . i.e. 



o = C 



1 + 



CO 



+ O(ln- 2 



(B.34) 



ln^/X 

where Co is universal and can be computed by using the perturbative expression of 
the scaling dimension of O. More precisely, if the scaling dimension 7e>(w,f) has the 
perturbative expansion 

lo(u, v) = g 00 + go\v + g 2V 2 + g 2 ou 2 + . . . (B.35) 

we obtain 



co 



-ngo2- 



(B.36) 



Corrections proportional to 1/ In £/X should instead be absent in RG invariant 
quantities. Indeed, if R is such a quantity, if we neglect the scaling corrections, R 
satisfies the scaling relation 



R(t)=G r [Q,v( 



(B.37) 
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for any I. This implies that R(t) is independent of v(l), hence an analytic function of 
Q and therefore of r. It follows 

R(r) =R* + + 0(ln" 4 £), (B.38) 

In i/X 

where the costant cr is expected to be universal. 

Appendix C. General behavior close to a critical point 

Let us consider a multicritical point in a two-parameter space labelled by T and a and 
let us assume that the correlation length behaves as 

£(T,a)~[T-T c (0)]^ a = 0, (C.l) 

£(T, a) ~ [T - T c {a)\- V * a > 0, (C.2) 

where T c (cr) is the u-dependent critical point and V\ ^ u 2 . According to the RG, close 
to the multicritical point £ (T, a) behaves as 

C(T,a) = u t (T,a)-^F[u a (T,a)u t (T,a)-% (C.3) 

where u a (T,a) and u t (T,a) are the scaling fields and <fi and u m two critical exponents. 
Since one of the two scaling fields must vanish along the transition line, we define u t (T, a) 
as the scaling field which has this property. Therefore, we define 

-( T -") = Z W i - (a4) 

For a — > and T — > T c (0), it behaves as 

u t (T, cr) = r + c a a + . . . r = T - . (C.5) 

We assume that c CT 7^ 0, i.e. that the transition line is not perpendicular to the line 
a = 0, as it occurs in the RPXY model. Finally, we note that u a (T, a) does not vanish 
on the transition line, unless a = 0. 

Now consider T — > T c (cr) at fixed novanishing a. Since u a (T,o~) ^ we obtain 
flCj) only if 

F{x) ~x x A = (C.6) 

for x —7- 00. To go further let us distinguish two cases: (i) u a (T,a) vanishes identically 
for cr = 0, i.e. u a (T, 0) = for any T; (ii) u a (T, 0) is different from zero unless 
T = T c (a = 0). 

In case (i) (1C.1|) requires 

F(0) ^ 0, z/ m = v x . (C.7) 

Assuming ?v(T = T c (0), a) = rf^cr for cr — >■ we obtain 

i{T = T c (0), a) = (c a a)-^F(d a c-f'a 1 -f'). (C.8) 
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The observed behavior depends on the value of 0. For < 1, since F(0) 7^ we obtain 

£(T = T c (0), a) = (c a a)-^ (a + 6a 1 -* + . . .) (C.9) 

The corrections are correct provided that F(x) is analytic for x = 0. If > 1, using 
(1C.6I) we obtain the behavior 

£(r = T c (0), a) ~a- v V = v x - (1 - 0)A = ^ ~ X ) + ^ . (C .l ) 

In case (ii), if w CT (T, er = 0) = g^t + 0(r 2 ) we obtain for a = 

f(T, 0) = T- Um F(d T T l - (t> ), (C.ll) 
which shows that 

F{d T T l -' t> ) ~ r 1 ^ 1 (C.12) 

in the limit r — >• 0. Let us now consider the behavior for T = T c (0) as a function of a . 
For a — > we have 

£(T = T c (0) )( x) = c a Vm a- Vm F{d a c-J'a 1 ^) ~ a^ 1 , (C.13) 

where we have used relation (IC.12I) . Thus, in case (ii) we have £(T = T c (0),cr) ~ o"-^ 1 
for any value of 0. 

Let us now show that the case relevant for the RPXY model is case (i). Indeed, 
case (ii) can only occur if the two relevant operators which occur at the multicritical 
point are both present in the model at o = 0. This does certainly not occur in our 
case in which a is associated with randomness. Therefore, our result that in the RPXY 
model £(T = T c (0), a) behaves as a~ ui implies that <fi < 1, i.e. that the RG dimension of 
the new operator that arises in the theory with a ^ is less relevant than the thermal 
operator present at a = 0. This is also the case of three-dimensional randomly dilute 
Ising systems or ±J Ising models at their ferromagnetic transitions at small disorder. 
Indeed, the crossover from the pure critical behavior to that of the randomly-dilute Ising 
universality class is described by the crossover exponent = a>i s = 0.1096(5) [69, 80J, 
see also the discussion reported in [SI] . 

Similar considerations apply to other quantities. For instance, consider a RG 
invariant quantity R. It behaves as 

R(T,a) = r[u (T {T,o-)u t (T,a)-' l> \. (C.14) 

If < 1, R(T,o~) approaches the same value R* along the lines o = and T = T c (0). 
Morover, in the second case we expect corrections of the form 

R(T e {0), a) = R* + aa 1 ^ + . . . = R* + a'^~ 1)/ui + ... (C.15) 
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Appendix D. RG equations in the presence of randomness 

The RG equations in the small disorder regime and close to the paramagnetic-QLRO 
transition line have been derived in [3j [231 1211 EU [35] : 

dT . 3 , r2 

a = ' 

da 
dY 

_ = ( 2 - tt/3 + naP 2 )Y, 

where Y is the vorticity and only terms up to 0(Y 2 ) are kept. Let us now redefine the 
couplings as follows: 

T- 1 = -{2 + v + a), 

For u, v — y the RG equations become 
du 

a = " W) 

dv 2 

da ,„ . 

Ti=°- < D ' 2) 

We have thus reobtained the RG equations for the XY model. This implies that, in the 
region of couplings in which ( ID. 21) hold, the RG behavior is analogous to that close to 
the KT fixed point, apart from an analytic redefinition of the scaling fields. 



Appendix E. Magnetic correlations in the gauge-glass model 

For the gauge-glass model (a = oo) we can derive some identities which relate magnetic 
and overlap quantities. The basic observation is that for a = +oo the distribution 
function of the A xy variables is gauge-invariant. Hence we have 

i(r xi . . . ca> ■ ■ ■ vw>] = k ■ ■ ■ v iv yi ■ ■ ■ v yn [« . . . r x j yi ■ ■ ■ vosu) 

for any set of phases V x . It implies that magnetic correlations vanish unless each Xi is 
equal to some yj. Analogously we have 

[« . . . r x Jyi ■ ■ • ■ ■ ■ <A • • • M] = (E.2) 

= K ■ ■ ■ v L v v, ■ ■ ■ v yn v z : . . . v;v tl ...v tn [{r xi . . . r x J yi ■ ■ ■ Vw> « • • • • • • VOl- 

These relations allow us to write 

WM] = *xv, (E.3) 

\( 1 Px 1 } Px2 l i J yi ~~ ^xiyi^X2V2 $xry<2$X2yi ~ $xryi$xxX2$xxy<2i (E-4) 

[(^x^yi) (V^^Wl = fix 1 y 1 $x 2 y2 + $Bijn$Bayi [| (VVi^Vi) I ] — fix 1 y 1 fix 1 x 2 fixl$-5) 
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It follows 



[<H 2 >] = v, 

[(\^}} = 2V 2 -V, 

[M 2 ) 2 } = V 2 + V\ -V, (E.6) 



which imply 



X = l, 

X4 = 1 - 2xo, 

X22 = Xo ~ 1- (E.7) 

Moreover, it is easy to show that £ = 0. Relations ( IE. 71) show that %4 an d X22 both 



diverge as Xo- I n the critical limit we have Xo ~ £0 because ?7 = 0. Therefore we can 
write 

X 4 ~ ~2ail X22 ~ <o, (E.8) 

for £ a — )• 00, where a is constant. 

We shall now assume that these results are valid for the whole universality class: 
for any o > ap, relations (1E.8j) always hold with a constant a which in general depends 
on a. We can reexpress these results in terms of the quartic couplings. If we use (IE. 81) 
we have 

x 2 e 
x 2 e 

Since the magnetic susceptibility x an d correlation length £ are finite and nonzero 
(except for o = 00, where anyhow the quartic couplings are not well-defined since 
£ = for any L), we expect that and (722 diverge as ^ in the critical limit. As for 
g c = #4 + 3(?22, ( IE. 101) shows that the leading ^ term cancels. Since in the calculation 
we have neglected the scaling corrections to ( IE. 81) . this does not necessarily imply that 
g c remains finite in the critical limit, but only that g c £~ 2 — > as £ D — > 00. The exact 
behavior depends on the neglected scaling corrections. These predictions are confirmed 
by our numerical results, see Sec. 16.71 It is worth mentioning that this behavior is 
analogous to that observed in the 2D Ising spin glass model, where %4 behaves as Xo 
and thus diverges approaching the glassy transition; see, e.g., [82J and references therein. 



922 = (E.10) 
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